Neutral test particle orbits in the Kerr— Newman spacetime 



Daniela Pugliese^, Hernando Quevedo^'^'^, and Remo Ruffini^ 
^ School of Mathematical Sciences, 
Queen Mary, University of London, 
Mile End Road, London El ^NS, United Kingdom 
^ Dipartimento di Fisica and ICR A, 
Universitd di Roma "La Sapienza", 
Piazzale Aldo Moro 5, LOO 185 Roma, Raly 
ICRANet, Piazzale della Repubblica 10, 
L 65 122 Pescara, Raly 
^Instituto de Ciencias Nucleares, 
Universidad Nacional Autonoma de Mexico, 
AP 70543, Mexico, DF 04510, Mexico 
"^Instituto de Cosmologia, 
Relatividade e Astrofisica and ICRANet - CBPF 
Rua Dr. Xavier Sigaud, 150, 
CEP 22290-180, Rio de Janeiro, Brazil 
(Dated: March 26, 2013) 

We present a detailed analysis of the orbital circular motion of electrically neutral test particles 
on the equatorial plane of the Kerr-Newman spacetime. Many details of the motion in the cases of 
black hole and naked singularity sources are pointed out. We identify four different types of orbital 
regions, which depend on the properties of the orbital angular momentum, and define four different 
kinds of naked singularities, according to the values of the charge-to-mass ratio of the source. It is 
shown that the presence of a particular type of counter-rotating test particles is sufficient to uniquely 
identify naked singularities. It is pointed out that the structure of the stability regions can be used 
to differentiate between black holes and naked singularities. 
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I. INTRODUCTION 



The Kerr-Newman (KN) spacetime is an exact solution of the Einstein-Maxwell equations that describes the exterior 
gravitational and electromagnetic fields of a rotating charged source with mass M, angular momentum J and electric 
charge Q. In Boyer-Lindquist coordinates, the KN line element can be written as 



( A - sin^ (9) 2asin^(9(r2 



A) 



(r^ ^a^y - Aa^ sin"' 6> 



where we used geometrized units with G 



and 



1 = c, and 



sm 



-dr' 



dtd<p+ 



2Mr 



+ cos^ 9. 
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The parameter a stands for the angular momentum per unit mass, as measured by a distant observer. The corre- 
sponding electromagnetic vector potential 



Qr 



[{dt)a - asiii^ d{d(j))a] , 
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depends on the charge Q and the specific angular momentum a. It then fohows that the magnetic field is generated 
by the rotation of the charge distribution. 

The limiting cases of the KN metric are the Kerr metric for Q = 0, the Schwarzschild metric which is recovered 
for a = Q = 0, the Reissner-Nordstrom (RN) spacetime for a = 0, and the Minkowski metric of special relativity for 
a = Q = M = 0. The KN spacetime is asymptotically fiat and free of curvature singularities outside a region situated 
very close to the origin of coordinates. 

Several critical points characterize the geometric and physical properties of this spacetime. In particular, the 
function A vanishes at the radii 



r± = M± - a2 - Q2 (5) 

which are real only if the condition > -\- is satisfied. In this case, r+ and r_ represent the radius of the 
outer and inner horizon, respectively, and the KN solution is interpreted as describing the exterior field of a rotating 
charged black hole (BH). In the case < + Q^, no zeros of A exist and the gravitational field corresponds to 
that of a naked ring singularity situated at 

+ cos^ = 0. (6) 
The norm of the timelike Killing vector field = for the KN metric in Boyer-Lindquist coordinates ([T]) reads 

A-a^ sin^ 

9tt = ^ , (7) 

Then, the norm is positive in the region where 

+ cos^ (9 + - 2Mr < 0, (8) 



with rf = M± - Q2 _ ^2 cos2 (9) 

where — — o? cos^ ^ > is satisfied. In particular, for a black hole it is r+ < r < a region which is known 
as the ergosphere. In the regions where becomes spacelike, the coordinate t cannot be used as a time coordinate 
and the analysis of the physical properties of this spacetime requires a different approach. 
The extrema of gu are located in 



* _ ± V2a2M2 + Q4 ^ 2a2M2 cos(2^) 
- 2M ' ^^^^ 

which is a function of the parameters a and Q, and the coordinate 9. At these radii we find that 



2a2 + (q2 -p ^2o2M2 + Q4 + 2a?M'^ cos(26»)) sec^ 9 



- 9u{rl) = ^ , (11) 

SO that Qtt has a maximum at r = r^. The limiting cases of these radii are of interest. When a = (Reissner-Nordstrom 
spacetime) we have 

rX = j^. rl=0, gu{rl) = -l^^. (12) 



For a Kerr spacetime, Q = 0, we have that 

= V cos2 ^, = — V cos^ ^, guir"^) ~ ~ 



(13) 



Va2M2 cos2 



In this work, we will perform a detailed analysis of the motion of test particles along circular orbits on the equatorial 
plane of the KN spacetime. Since test particles moving along circular orbits are particularly appropriate to measure 
the effects generated by naked singularities, we will focus on the study of the physical differences between a black hole 
and a naked singularity see also [THS] and for further applications [H |5] . 

This work follows a series of papers 0-[8J in which we showed the presence of a typical band structure, i.e., 
disconnected regions of stable orbits for the spacetimes generated by Kerr and RN naked singularities. This structure 
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is absent in the case of black holes. It was also found that there exist two types of singularities that affect the 
properties of the orbits and their stability. We generalize here those works considering a field source with charge and 
intrinsic angular momentum. 

The plan of this paper is the following: In Sec. [llj we introduce the effective potential for a neutral test particle 
moving along a circular orbit on the equatorial plane, and we set the main notations to be used in this work. In 
Sec. |III[ we explore the motion in the spacetime of a KN black hole while the case of a naked singularity is explored 
in Sec. |IV| Section [V| summarizes the main results and compare the two cases. Finally, the conclusions in Sec. |VI| 
close this article. 



II. CIRCULAR ORBITS 

We consider the circular motion of a test particle of mass /i in the background represented by the KN metric ([T]). 
We limit ourselves to the case of orbits situated on the equatorial plane. We adopt the effective potential approach 
to the study of the test particle dynamics. The test particle motion is therefore described as the one-dimensional 
motion of a classical particle in the effective potential V{r). In the case of equatorial geodesies {0 = 7r/2), one obtains 
the effective potential [9^ 

^ , - 4AC ^^^^ 

^ = -2A+ 2A - 

where 

A = (r2+a2)'-a^A, (15) 
B = -2aL (r^ + - A) , (16) 
C = a^L^ - {M'^r'^ + L^) A . (17) 

Here L = jigc^js^^^u^ is a constant of motion associated with the angular momentum of the test particle with mass /i 
and 4- velocity . The Killing vector = 9^ represents the axial symmetry of the rotating source. The additional 
Killing vector — dt is timelike and represents the stationarity of the field configuration. It also generates a constant 
of motion E = —iiQc^^^^^u^ which is associated with the total energy of the test particle, as measured by a static 
observer at infinity. 



Notice that the effective potential ( 14 ) was obtained originally as the solution of a quadratic algebraic equation so 
that the solution 



^ --2A 2A 

is also allowed. Nevertheless, this solution does not need to be analyzed separately, because its properties can be 
obtained from V by using the symmetry V~{—L) = —V{L). 

Circular orbits are given by the simultaneous solutions of the equations 

V = El II and ^ = 0. (19) 

dr 



The potential (14) is quadratic in the charge Q. We therefore will limit ourselves to the study of the case Q > 0. 
Moreover, the effective potential is invariant under a simultaneous change of the sign of a and L. Consequently, 
without lost of generality, we can restrict the analysis to the case a > 0, and will differentiate between corotating 
particles with L > and counter-rotating particles with L < 0. 

To perform a detailed analysis of the effective potential of a naked singularity with + > M^, we will consider 
the cases 

a > M, or > M, (20) 

otherwise 

M >a> as, and Q < M, = V^M^ - Q2 ^ (21) 

or 



M>Q>Qs, and a < M, = v^M^ - . (22) 



For a discussion of the motion in a KN spacetime see p^OHTS] and also [T9]-[2T]. 
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A. Sign of the effective potential V 

In this section, we explore negative energy states for circular orbits considering the sign of the potential given in 



Eq. 14 First, the effective potential is null {V = 0) for a particle with angular momentum L = Lq < with 




FIG. 1: The particle orbital angular momentum Lq/M/j. is plotted as a function of the charge-to-mass ratio Q/M and the intrinsic angular 
momentum a/M of the KN source for a = 0.71M (left plot) and a = M (right plot). 



Wm = -^\/-A/Ao, (23) 

where Aq = A(a = 0) = — 2Mr + Q^. This function is plotted in Fig. [l]in the regions showed in Tabl e [l} For a 
comparison with the most important limiting cases, we introduce also the radii = r±{a = 0) = M ± y^M^^^Q^ 
for the horizons of the Reissner-Nordstrom black hole, and = r±{Q = 0) = M ± \J — a? for the horizons of the 
Kerr black hole, respectively. Notice that coincides with the ergosphere boundaries rf on the equatorial plane 
= 7r/2. On this plane the ergosphere is completely independent from the source spin. In particular, for a KN naked 
singularity there is a change of sign in the time-time component of the metric tensor (and in the norm of the timelike 
Killing vector) inside the region [r^, r+], for a < M and Qg < Q < M and a > M with Q < M. 
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TABLE I: Table provides a description of the regions characterized by V/fi = with L = Lq. The particular cases of Kerr spacetime 
Q = 0, Schwarzschild (Sch) Q = a = 0, Reissner-Nordstrom (RN) a = and extreme Reissner-Nordstrom case (ERN) a = and 
Q = M are also listed. The outer r+ and inner r_ horizons, as well as r^^ = r± {a = 0) = M ± v'MS - Q2 for the horizons of the 
Reissner-Nordstrom black hole, and = r±(Q = 0) = M ^ V M'^ — for the horizons of the Kerr black hole are introduced. 



It can be shown that in general the effective potential is negative in the region L < Lq and r < r^j with 



4 4\/3M 4Q2 " 



(t1/3 73 



rl/3 



3ai/3 



a = (o^ + 18o (SM^ + 2Q'^) + Q\fiy'l&Q^ + (M2 + Q2) + ^2 (27M^ + 36M2Q2 + gQ^)) , 



/a4-2a2 (6Q2 + 0-1/3) + 0-2/3 



ri/3 



(24) 
(25) 
(26) 
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The behavior of this radius is depicted in Fig. |2] for different values of the intrinsic angular momentum and charge. 
It then follows that there exist different regions where the effective potential becomes negative; the results of this 
analysis are listed in Table [TH 



fV/M=constant rp-ZM-constant 




a/M a/M 

FIG. 2: The radius Va/M is plotted as a function of the charge-to- mass ratio Q/M and the intrinsic angular momentum of the Kerr- 
Newman source (left plot). The plane Q = M, a = M, and the outer and inner horizons are also plotted. The density regions ra/M = const 
are plotted for selected intervals of a/M and Q/M, in the region of a black hole (center plot) and in the case of a naked singularity (right 
plot). 
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TABLE II: Table provides a description of the regions characterized by V/fi < 0. The outer r+ and inner r_ horizons, as well as the 
radii = r±{a = 0) and = f^±iQ = 0) are introduced. 



Not all the regions in Table |n| are characterized by the existence of negative energy circular orbits. Indeed, this 
kind of orbits requires the fulfillment of the additional condition ^ = in those regions. Nevertheless, Table |ll| 
provides an indication of the regions where these orbits may exist. It is important to note here that, in general, 
negative energy states can exist for negative orbital angular momenta (i. e., counter rotating particles) and for 
sufficiently large angular momenta, that is, for L larger in magnitude than Lq. But the orbital regions in which this 
can happen strongly depend on the source charge, even in the case of neutral test particle orbits. Notice also that in 
the < Q < M case, the horizons for the Reissner-Nordstrom metric = rf" must be considered. Table [ll] shows 
the situation for the three cases Q = 0, a = 0, and (a 7^ 0, Q 7^ 0); the discriminant condition between the various 
cases is represented by the case Q > M, for which one has a naked singularity. In this case, the negative potential 
regions are located in a closed, bounded and continuous area of orbits between the singularity and the orbit r^r, the 
extension of this region will depend on a and Q as shown in Figs. [2] and, in general, it increases with both the charge 
and the spin of the source. 



6 



B. Circular motion around black holes and naked singularities 

In this section, we introduce the basic equations of the circular motion in the KN spacetime, we fix the main 
notations and introduce the most important quantities that are relevant for the characterization of the motion both 
in the black hole and in the naked singularity case. 

In order to explore the circular motion around black holes and naked singularities, we first study the condition 
y = 0. Thus, we analyze the equations for the circular motion 



F'(r,L,a,Q)=0, V = E/^, 



(27) 



and solve them with respect to the angular momentum of the test particle. We obtain that the general solution 
corresponds to L = where 



1 /s±2M2Vn 



with 



^ = 

n = 

s = 



(Q^ - rM) + [2Q^ + (r - 3M)r]^ 
-a" - rM) [a" + + (r - 2M)r]' [a" {Q^ - Mr) + {2Q^ - 3Mr) r^\ 



{-{Q^ - Mr) r^ [2Q^ + (r - 3M)r] + {Q^ - Mr) [2Q^ - (5M + r)r] + 
[2Q^ + Q^{r - llM)r - 2Q^(r - 2M)r^(5M + r) + 2Mr^[r(3M + r) - GM^]] } . 

The evaluation of the corresponding energies leads to 

(Q2 _ 2Mr) + Jr^ [a^ + . 



+ (r - 2M)r] (r^ [L\ + r^) + [r(2M + r) - Q^\) 



r^ + a2 [r(2M + r) 



(28) 

(29) 
(30) 

(31) 
(32) 



Now, the main point is to find the regions in the space of the radial coordinate and parameters {r /M^Q /M^a/M} 
where the solutions for the angular momentum and the energy are well defined. For instance, in Fig. [3j we show 
the behavior of the particle angular momentum and energy in the case a = I.IM for different values of the naked 
singularity charge Q. In Fig.[4j we present the numerical results for the case a = 2M and different values of the naked 
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FIG. 3: The angular momentum of circular orbits is plotted as a function of the orbital radius for different source charge-mass 
ratio Q/M. Circular orbits in the KN spacetime with a — I.IM are explored. Right: The energy V/ ii of circular orbits as a 
function of r/M for selected values of Q/M. 
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singularity charge Q. First, we study the orbital regions in which circular motion occurs, that is, where solutions of 
Eq. (27) for L exist. These regions have as boundaries the solutions of ^ = or explicitly 



(4Q2 + 9M2)r^ 



4 (3Q^ + a^) Mr + AQ^ + 4a^Q2 = 0. 



The solutions of Eq. ( 33 ) are the radii 

rt = - 9M 
^ 6 



/3T ± V3 W 27M^ - 24Q^ + - T 



(33) 



(34) 
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FIG. 4: Circular orbits in the KN spacetime with a = 2M are explored. The angular momentum of circular orbits is plotted as a function 
of the orbital radius for different source charge-mass ratii Q/M. Right: The energy Vf/j, of circular orbits as a function of r/M for selected 
values of Q/M. 
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where 



- 8Q 



A/3 



-1/3 



(8Q^ - 9M2) - 36a^ {27 M"^ - 42M'^Q'^ + 16Q 



(35) 



The study of the regions confined by these radii wih be the key point to analyze the circular orbits of neutral particles 
in a KN naked singularity and to distinguish this case from the black hole one. 

We see that in general there are four possible different radii and all of them depend on the value of the parameters 
M, Q and a that characterize the rotating source. For a given value of the source parameters not all the solutions 
lead to real and positive radii so that many different cases can arise. For a better presentation of the results it is 
convenient to introduce the notation 



ri = r_, r2=r^, ra = r^, r4 = r|. (36) 
The behavior of these radii is illustrated in Fig. ([5|. Notice that in general ri < r2 < ra < r4. The existence of the 




Q/M 



FIG. 5: Left: The radii ri (black surface), r2 (white surface), rs (gray, black-shaded, surface), and (black, white-shaded 
surface) plotted as functions of Q/M G (0,2) and a/M G (0,2); it is ri < r2 < rs < r4. Center: Surfaces (ri, r2, ra, r4) in the 
the region Q/M G {0, 1.1}, a/M G {0, 1}, r/M G {0,4}. The black region represents r < r+, where r+ is the outer horizon of 
the black hole. Right: The region Q/M G {1, 1.07}, a/M G {0, .251}, r/M G {1,2.35683} corresponds to the case of a naked 
singularity. 



surfaces is studied in detail in Fig. [6j The shaded regions contain all the points at which Eq. ( 33 ) does not possess 
real positive solutions. In the light gray regions of Fig. |6j four solutions exist, and in the white regions only certain 
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FIG. 6: Left: Plot of and as as a functions of Q/M. Right: Plot of Q4 (black curve), Qs (gray curve), Q5 (dashed cu rve), Qq (black 
thick curve) as functions of r/M. Black regions are forbidden: no solution of Eq. ( |33| exits. In the white regions, Eq. ( |33| ) is satisfied for 
r = rs and r = r4, and in the light gray regions, it is satisfied for r = ri , r = r2, r = rs and r = r4. In particular, for Q = and a = the 
solution is r = 3M, for < a < M it is r = r2 , r = rs and r = r4. For a = M it is r = M and r = 4M, finally for a > M it is r = rs = r4. 



radii, as indicated in the plots. Thus, to analyze these regions in detail, it is convenient to introduce the following 
charges 
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and the following spin parameter 



Q2 (8Q2 _ QM^y 



3V3M 



(37) 
(38) 

(39) 



which guarantees that in Eq. ( |34[ ) is real. 

Alternatively, it is also possible to investigate the solutions of the equation ^ = in terms of the source charge as 
a function of (a,r). In this way we will have a different view of the regions of existence of circular orbits. Thus, the 
solution of ^ = is 



(3M - r)M -a?—± ^Ja? {a? + 2(r - M)r) 



(40) 



The corresponding surfaces are studied in detail the left plot of Figs. [Tj which also includes the regions where these 
charges exist. It is convenient to introduce the following spin parameter 



and the radii 



= ^/2r{M-r), a^, = \ .J{r -ZMYr/M, 



^ = I (M±lVM^-2aA, 



(41) 



(42) 



which define the region r G]rM_,^M+[ where no solution of Eq. (33) exists (see the right panel of Fig. [7]-right). The 
radii 



TmJM 4 sin 



arccos (l - 2a^ /M'^^ 



JM 4cos 



arccos (2a 



r^jM = 2 1 + sin 



arcsin (l - 2a^/M'^^ 



(43) 
(44) 
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FIG. 7: Black regions are forbidden: no solution of Eq. \33\ exits. Left: Plot of an and as functions of r/M. The solutions are 
shown in the light-gray region and in the white regions. Right: Plot of (gray curve), (dashed black curve), (dotted-dashed 
black curve), rM_ (gray curve), tm^ (gray dashed curve) as functions of a/M. Black regions are forbidden: no solution of Eq. ( |33| > exists. 
In the white regions, the solution is Q = Q^; in the light gray regions the solutions are Q = 



were introduced in [7 , with a different notation, in order to study the existence of circular orbits in the Kerr spacetime. 
In Figs. [H 




FIG. 8: Left: Plot of (thick black curve), (thick gray curve), r* (black curve), re (dotted-dashed curve) as functions of Q/M. The 
solution of Eq. (|33) is in the white regions. Center: Plot of Qi and Q^jM — yJrjM as functions of the radius r/M. Dotted-dashed 
curve signs rg. In tne white regions, Eq. ( |33| is satisfied for a = ai. The plot is therefore divided in the regions: < r < M, M < r < 3M, 
r > 3M. In particular, at r = M and Q = 1, a solution exists for a > 0. The dotted-dashed line is Q = M, the dashed line is r = M and 
the dotted line is r = 3M. Right: The black hole angular momentum ai is plotted as a function of r/M, and the charge-to-mass ratio 
Q/M of the black hole. The outer horizon plane r+ (gray) is also plotted. The black region is r < r+. The curve r = (black thick) is 
also plotted. 



we show the solution of the equation ^ = in terms of the spin parameter 



where we used the definitions 
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1 [2Q^ ^ {r - 3M)rY .^^ 



and Qi/M = \f^^^^ with Q^/M = ^Jr/M (that corresponds to the definition = Q'^ /M). The radii (r^, r*), 
introduced in [6], characterize the dynamics around a Reissner-Nordstrom naked singularity {Q/M > 1) where the 
following inequality holds < r~ < r+, and are the radii at which the value of the angular momentum and the 
energy of the test particle diverge. At the classical radius circular orbits exist with "zero" angular momentum. 
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This means that a static observer situated at infinity would interpret this situation as a test particle that remains 
motionless as time passes. This phenomena can take place only in the case of a naked singularity and is interpreted 
as a consequence of the "repulsive" force generated by the charge distribution. It is interesting to note that these 
radii have an important role even in the presence of a non zero source spin. In particular, it turns out that from 
Eq. (28) the orbital angular momentum L is defined when 11 > 0, and from Eq. ([30|) we recover r > r>^, at r = it is 



l^^\J Q'^{Q'^-M'^) -> which is in fact zero in the Reissner-Nordstrom case. In Fig. |8[ we plot the behavior of these 
radii as functions of the ratio Q/M and a/M. 



Now, for the existence of solutions of the Eq. (27) it is however necessary to demand that the solution (28) be real. 
Therefore, we study the solutions of the equation T? — AM^Ii = 0, or 



Q^a^ + {4Q^a^ + a^) Mr - 2 {2M^ + Q^) aV^ + 2MaV^ - QV^ + Mr^ = 0, 



(47) 



which will be denoted by r (see Fig. [9| . But in general we can also find the regions of circular motion by fixing the 
value of the source angular momentum a and plotting the ratio Q/M in terms of the radial coordinate. Thus, we can 
also express the solution of Eq. (47) in terms of the charge Qc 



Qc 



\ 



-a^ + 0? (-2(r - 2M)r + a'^ + 2a? {r - 2M)r + r^^ + r"^ i^-r"^ + a'^ + 2a? {r - 2M)r + r"^^ 



2a? 



• (48) 



In this case, it is convenient to introduce the charge 



Qr = V-a2 - (r - 2M)r , (49) 
which is the solution of the equation A = 0, and the radius Tx which represents the real solutions of Eq. (HtI when 
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FIG. 9: Left: Plots of r as functions of Q/M and a/M. The black region is inside the horizon r < r+. Center: Plot of tq in the case 
Q = and Q = M (inset plot). Right: Plot of in the range [0, 3M) as a function of a/M. The inset plot shows a zoom for a < M. 



the source charge is Q^^ or 



12Ma'^ + (4a^ - M^) r - 12MaV^ + 5aV^ - 2Mr^ + 2r^ = . 



(50) 



This radius is plotted in Fig. [9] for different values of the ratio a/M. Notice in particular that Tx does not depend on 
the source charge but only on its angular momentum. 
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III. BLACK HOLES 



In this section, we study the Kerr-Newman black hole case, considering the constraint (5 > 0, where S = M^—Q^—a^. 
In particular it is (5 = for + = M^, and Q G [0, M] and a G [0, M]. For Q = a = 0, Schwarzschild black hole, 
it is (5 = M. 




FIG. 10: The outer horizon r+ = M + — — Q^, the outer ergosphere radius on the equatorial plane r'^ = M -\- sj — 

and the radius rj" = 2r^, are plotted as functions of the black hole charge-to-mass ratio Q/M, for different values of the intrinsic angular 
momentum a/M. 



In Fig. [Toj the outer r+ and inner r_ horizon surfaces are plotted as function oi a/M and Q/M. We are interested 
in the solutions of Eq. (27) in the region r > r+. Solving Eq. ( [27| ) for the angular momentum L, we obtain the 
solutions given in Eq. (28). In Figs. 11, the angular momentum L/^Mfi) of circular orbits is plotted as a function of 



the orbital radius for different charge-to-mass ratio of the black hole and its intrinsic angular momentum. There are 





FIG. 11: The angular momentum L/(m/i) of circular orbits is plotted as a function of the orbital radius for fixed source charge-mass 
ratio Q/M and intrinsic angular momentum a/M^ with M < a < \/l — Q^. 
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zones in which two kinds of orbits, characterized by two different angular momentum, but same orbital radius, are 
possible. In particular, it is possible to show that there are no circular orbits in the KN black hole spacetime with 
zero angular momentum, i. e., there are no solutions of the following equations 



V\r,Q,a,L) = 0, L = 0. 



(51) 



However, one can see that V {r^Q^a^L = 0) > in the region r > r+; this means that for particle with an "angular 
momentum" L = the effective potential monotonically increases as the outer regions of spacetime are approached. 



see Fig. 12 



:0.9M 




FIG. 12: Left: The effective potential V/ for a Kerr-Newman black hole is plotted as a function of r/M and a/M, for different values 
of Q/M and orbital angular momentum L = 0. The lower plane (black) corresponds to Q = O.OIM, the middle one (white) to Q = 0.5M, 
the outer one (gray) to Q = M. The function V monotonically increases with r/M. The black hole angular momentum ai is plotted as 
a function of r/M, for different values of the charge-to-mass ratio. Center-left: The outer horizon r+ = M + ^/ M^ — o? — ^ the outer 
ergosphere radius on the equatorial plane = M + \J M^ — and the radius = 2r^ are plotted as functions oi r/M and the black 
hole charge-to-mass ratio Q/M. The angular momentum as = \J M^ — is plotted (gray curve); for the black hole case it is a < as- 
Left panel: Q = 0.5M , as = 0.866M, r+ = 1.866M, r+ = 3.732M, r~ = 0.268M. Right panel: Q = 0.9M, as = 0.436M, r+ = 1.436M, 
rj" = 2.872M, = 1.1282M. Shaded regions are forbidden: the horizon constraint of the KN black hole implies < a < as < 1 and 
r > r+. 



But before we deal with the general case a 7^ and Q ^ 0, it is useful to summarize briefly the known outcomes 
for the two limiting cases of Kerr {Q = 0) and RN (a = 0) spacetimes. We refer for details to the extensive literature 
on the subject and to the aforementioned works [BHHl- In fact, we will see that the source spin and a charge deforms 
the regions of circular motion and the stability properties of a RN and Kerr spacetime, respectively. To discuss the 
dynamics of a black hole spacetime it is convenient to introduce here the following radius 



2(M± -Q2)^ 

We notice that = 2r^ and, for circular motion on the equatorial plane {6 = 7r/2), it is 



Figs. |1Q|12| Furthermore, it is possible to prove that the following limits are satisfied: 

= ^3, for (5 = (Kerr metric), 

= r3 = r4, for a = (Reissner-Nordstrom metric). 

Moreover, following the notation of |7], we define the radius = r4((5 = 0), i. e.. 



2^RN 



2r? 



(52) 



(53) 
(54) 



Ta = 4M COS 



6 



arccos 



(55) 



Tables III and IV show the case Q = from two different points of view. This differentiation in the presentation 
of the orbital regions is important for the study of the stability problem also in the general case of a KN spacetime. 



In Table III, we consider different values for the source angular momentum a, ranging from the Schwarzschild black 
hole, a = 0, to the extreme Kerr black hole with a = M. Each range is then divided in various spatial regions, each 
of which is characterized by a particular value of the orbital angular momentum L of the test particle. 



In Table IV, we present an alternative classification of the circular orbits. The radial coordinate is split into various 
regions that range from r = M to r ^ 00. It turns out that the radii r = 3M and r = 4M are of particular interest 
because they determine spatial boundaries in which the value of the orbital angular momentum L drastically depends 
on the source angular momentum a. In each spatial region, we present the type of allowed circular orbits (rotating 
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Case: Q = (Kerr) 



a = 


(Schwarzschild) 


<a < M (Kerr) 


a = M (extreme Kerr) 


Region 


L 


Region L 


Region L 


r > 3M 


±L- 


]rM3,ra] L- 




]M,4M] L_ 






]ra, oo[ {L-,-L^ 


-) 


]4M,oo[ {L-,-L+) 



TABLE III: Angular momentum of circular orbits L in the region r > M of a Kerr spacetime for specific values of the source angular 
momentum a. 



Case: Q = (Kerr) 



M < r < 3M 


3M < r < 4M 


r = 4M 


r > AM 


a 


L 


a 


L 




a 


L 


a 


L 


]ai,M] 


L_ 


[0,a4 
[a^,M] 


(i>- 
L_ 


-i>+) 


[0,M[ 
M 


(i^-,-i^+) 
13/{4.V2)Mfi 


[0,M] 


{L-,-L+) 



TABLE IV: Angular momentum L of circular orbits in the region r > M of a Kerr spacetime for specific values of the radial coordinate 
r. 



and counter-rotating) with the corresponding value of the angular momentum of the test particle {L±). More details 
about this classification can be found in [6l [7j • Finally, in the limiting case of a RN black hole (a = 0) , one can see 
that circular orbits with L = zbL_ are allowed only in the region outside the radius r+. 

The main question now is to determine how the electric charge of the gravitational source will affect this classifi- 
cation. Since the charge generates its own gravitational field, it is expected that its presence will modify the spatial 
regions as well as the ranges of the source angular momentum in which circular motion can take place. 
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A. Orbital regions around a KN black hole 



We give here a description of circular orbits with L 7^ around a KN black hole from the outer horizon r+ to 
infinity. Since it is always r+ > M, in order to simplify the exploration of the dynamics in the region r > r+, we 
approach the more general problem of finding the extremes of the effective potential in the larger region r > M with 
the constraints < a < and < Q < M. The results of this analysis are summarized in Tables [V| and |VI[ Notice 
that the spin parameter a = plays an important role for the definition of the different ranges of Table. [Vj because 
for Q > M a black hole case occurs if the source spin satisfies the condition a < a^. We point out that there are two 
different orbital regions, namely, the region I, with L = L_, and the region II, with L = Table |V| shows 

the allowed angular momenta of the test particle (L), depending on the the values of source spin (a), for fixed ranges 
of the charge Q/M and different orbital regions defined by r/M. We also include in this table the particular case of 
an extreme RN black hole, Q = M and a = 0. This is a generalization of the Table [TV| for the case Q ^ 0. We can 
see t hat the regions in which circular motion is allowed do not differ qualitatively from the configuration shown in 



Tab. 



IV 



for the case Q = in which the limits 



r+ - 



3M and 



4M hold. We conclude that the introduction of 



a source charge does not change qualitatively the structure of the orbital regions, but only their extensions, i.e., the 
boundaries become modified as shown in Table \V\ 

The extensions of the orbital regions I and II will depend depend ona/M and Q/M but, at least for the characteristic 
values of black hole spacetimes, we do not see the appearance of an additional type of orbit. For example, there are 
no counter-rotating orbits with L = — L_ in r > r+ as in the limiting case a = 0. 



Case: < Q < VSM/2 {V3M/2 < Q < M) 




{r^)M <r <r^ 


< r < 


r = r^ 


r > r+ 


a L 


a L 


a L 


a L 


]ai,as\ L- 


[0,ai[ (L_,-L+) 
[ai,as] L- 


[0,as[ {L-,-L+) 
as L- 


[0,as] {L-,-L+) 


Case: Q = M 




r>2M 


a L 


a = ±L_ 



TABLE V: Description of the circular orbits for a test particle in a Kerr-Newman black hole with < Q < M, in the region r > M. The 
orbital angular momentum is given for each region. This is a generalization of Table [Tv] to the case Q ^ 0. 



The results presented in Table [V| are illustrated in Figs. [T3| and [14 

Case: < Q < M 



a = 


(Reissner Norstrom) 


< a < as 


a = 


= as 




Region 


L 


Region 


L 




Re^ 


lion 


L 


r > 




]r3,r4] 


L- 




]r3, 


r4] 


L- 




]r4,oo[ 


(L- 




]^4, 


oo[ 





Case: Q = M 



a = 



(Reissner Nordstrom) 



Region 



L 



r>2M 



TABLE VL Description of the circular orbits for a test particle in a Kerr-Newman black hole with < Q < M, in the region r > M. 
The orbital angular momentum is given for different values of the black hole spin a/M. This is a generalization of Table [HII to the case 
Q/0. 



On the other hand. Table VI shows, for different intervals of the specific charge Q/M, the regions of existence 
of circular orbits for different ranges of the source spin a. For this analysis, it is convenient to define the following 
angular momentum parameters 



/ _ 2r(r - 2M)(rM T Q^) + (2rM T Q'^)\/Q^ + 4Q'^M(r - M)r - 4(r - M)Mr'^ , , 

" V 2(rM-Q^) ' 

which determine the boundaries of the different regions. In the limiting case a+ = a_, these regions reduce and very 
particular orbits appear which are studied in detail in Appendix [A] 
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a=0.1M 



a=0.8M 




a=0.9M 



(BH) 



i 



2 
r/M 




LJ-L_ 



^ 0.6 
) 

0.5 
0.4 

0-2 
0.1 

J 0.0 
4 u 

a=0.99M 





(BH) 














i . 


; (x+,-i-) 



2 3 
r/M 



(-i+i-) 

(BH) 



r/M 

FIG. 13: Black hole case: circular orbital regions in the plane (Q/M, r/M) for fixed values of a/M. The angular momentum of the 
orbiting test particle is shown in each spatial region. The curves are also plotted. See also Table V The horizons r± are 

represented as white curves. The dotted-dashed curve is the outer ergosphere rt • 



The orbital region boundaries are determined by the radii r+, rs and r4, which for fixed Q/M are in general 

Q=(V3"/2) M 



functions of a/M. The results of Table VI are illustrated in Figs. \\A 
Q=0.1M 



Q=0.5M 




FIG. 14: Black hole case: circular orbital regions in the plane {a/M^r/M) for fixed Q/M. The angular momentum of the 
orbiting test particle is shown in each spatial region. The intrinsic KN angular momentum a-/M (gray curve), a+/M (black 
curve) and aj / M (dashed gray curve) are also plotted as functions of r/M. The horizons r± are represented as white curves. 



See also Table 



VI 



The dotted-dashed line is the outer ergosphere . 



Finally, we notice that in Figs. 13|14 we have represented the entire orbital range r > 0, although for the study of 
circular orbits we are restricting ourselves to the physical region outside the outer horizon r > r+. Nevertheless, it is 
interesting to mention that, according to Figs. 13|14[ there exist circular orbits with L = — L_ in the region r < V-. 
The study of this motion inside the inner horizon is outside the scope of the present work. 
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B. Stability of circular orbits 



The radius of the last stable circular orbit for a neutral particle orbiting on the equatorial plane of a KN black hole 
is defined as the solution of the following equations: 



E, F'(r;,eo) = 0, V"{ri.,co)=0, with 



(57) 



as 



Stable (unstable) orbits are characterized by the condition V"{risco) > {V^\risco) < 0)- In Figs. 15, we plot ri 
a function of the intrinsic angular momentum a and the charge Q. We note that the radius risco is given as a surface 
that has as a limiting case the planes identified by the radii rff^^^ the last stable orbit for a = 0, and rf^^, the last 
stable orbit for Q = 0. The exact expressions for these radii were given in [7^ and [6], respectively. We can see from 



Figs. 16 that the surface contained between the two limiting cases is reduced to a curve if the condition 



is satished. This happens when Q and a are linked by the following relation 

r^- r 

9M 



RN 

Isco 



8M 



15M 



(58) 



In particular, if a = 0, we obtain that 4M < risco ^ 6M and if < a < M, then the last stable orbit in the spacetime 
of a KN black hole is located within the bounded region r/3 < risco ^ 6M where 

/ r I 1 

3a2 



rp/M 



+ a^M^ (5M2 + 2a?) + ^J a? {a^ - M^f (M2+4o2) 

l73~ 



M6 + a2M2 (5M2 + 2a2) + ^a^ {a^ - M^ f {M^ + 4a^) 



(59) 



This behavior is depicted in Fig. 16 It is interesting to mention the identity r/3 = r-^ isco 

at the boundary the effective intrinsic angular momentum ^Jl — [a/MY plays the role of the specific charge Q/M. 
Moreover, for the particular value a = 0.430075M we obtain that rf3 = r 



See Figs. 
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(Q = ttg), which means that 
3ific charge Q/ 
and|l6l Finally, for a = 





FIG. 15: Left: The outer horizon (black surface), the surface risco(L-) (hght-gray, meshed surface) and risco( 
black thick mash). Note that risco{~L-^) < risco{L—). Center: The outer horizon (black surface), t^^^ for a 
for Q = (white surface) and r^^^ for Q = (upper black-thick-mashed surface). Note that r^^^ = r^~^ at r = 6M and 
along the curve Q = Qk (see text). Right: The outer horizon (black surface), rf^^^ (gray surface), r^~^ (light-gray surface), / ^^^^ 
black-thick-mashed surface), the surface risco{L—) (light-gray, meshed surface) and nsco(~^+) (gray surface with black thick mash). 

and Q = we have that risco = 6M. 

The energies of the circular orbits are studied in Figs. 



L+) (gray surface with 
: (gray surface), r^^~ 

K- ^ ^RN 
Isco Isco 

(light-gray surface), r,^+ (white- 



17 In particular, we note that E[—Lj^) > E{L_) in the 



regions where both of them are allowed. As a function of the radial coordinate, the energy of circular orbits shows 
points of minimum values, i.e., stable orbits. There are also radii (rs and r^) at which the energy necessary for a test 
particle to remain in orbit around the source becomes infinite. This kind of circular orbits is also found in the RN 
and Kerr spacetimes. This fact is confirmed in Figs. [18] where, following the orbital classification of Sec. [Ill A[ we 
give a classification of the orbits energy for all radial regions. For some particular orbits (r = 4M, r+ and rj), the 
behavior of the energies is shown explicitly. 



The stability of these orbits can be studied by the comparison with the Kerr black hole case. In Figs. 19, the results 
are shown for black holes with different charges. Following [7 , we define two different radii for the last stable circular 
orbits as follows 



Isco 



risco{-L^), r 



Isco 



risco{L-) 



(60) 



r,(Q=0) 
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a/M a/M 



FIG. 16: Left: Last stable circular orbit radius risco for a Kerr black hole (Q = 0) - rj^"^^ (gray curve) and r^^^^ (black curve) -, the 
outer horizon r+, and the radius (dotted-dashed curve) are plotted as functions of the KN black hole spin a/M. The dashed line 
represents a = 0.430075M. Right: The charge Qk (black curve) of a KN black hole is plotted as a function of the momentum a/M, where 
^Isco ~ ^Isco (^^^ text). A naked singularity exists in the region Q > Qs- 



Q=1/2M Q=V3"/2M Q=9/10M 




r/M r/M r/M 

FIG. 17: Energy Ef/j, of circular orbits in terms of r/M for the orbits with angular momentum — and L_, respectively, for fixed Q/M 
and different a/M. The dotted line is E = fi. The radius rs (dotted-dashed line) and r4 (dashed line) are also plotted. The points on the 
curves represent the minima of these energies. 



for the orbits with angular momentum L = — and L = respectively. We now compare the general configuration 



for different values of Q/M ^ with the configuration in the case of a Kerr black hole ((5 = 0). We show in Fig. 19 
that for a KN black hole the orbital regions of stability are qualitatively similar to those of a Kerr black hole. The 
main feature is that if the charge-to- mass ratio increases, for fixed black hole spin a/M, the radius of the last stable 
circular orbit r^^^^ decreases. In particular Figs 20 show that, for fixed a/M and for different values of Q/M^ the 



stability regions bounded by risco{—L-^) and risco{L_) are connected regions. This is the main result of this section 
and it will be important when comparing with the case of naked singularities. The right plot of Figs. 20 clarifies the 
arrangement of the stability regions in a schematic way: in the region r+ < r < risco{L_) all orbits are unstable, in 
the region risco{L-) < r < risco{—Lj^) the orbits with L_ are stable while the orbits with L = — are unstable, 
and finally in the region r > risco{—L-^) all orbits are stable. This scheme does not change significantly when the 
parameters Q/M and a/M are changed; however, we note that the radius risco{L-) can lie on the right of the region 
I so that all the orbits with L = L_ are unstable. This occurs for values of a/M sufficiently small. As a/M increases, 
the radius risco{L-) becomes smaller and approaches the horizon, reducing in this way the regions of instability for 
the orbits L = L_ which belong to both the I and II regions (see Fig. |20]left). Furthermore, the radius risco{L_) 
splits the I into two different stability zones: Orbits with r > risco{L-) are stable whereas orbits with r < risco{L-) 
are instable (see Fig. [2Q| center). 
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Q=0 Q=0.5M Q= V3"/2 M 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 0.0 0.1 0.2 0.3 0.4 0.5 

a/M a/M a/M 



FIG. 19: The radius t^^^^ = risco{~Lj^) and t^^^^ = fisco{L—) of the last stable circular orbits for a KN black hole, with different fixed 
charge-to-mass radio Q/M, as a function of the black hole spin < a < a^. The horizon r+, the radius rZjr (dotted line), rj" (dashed line), 
ra (black curve), r4 (gray curve) are also plotted. The black region (r < r+) is forbidden. The value of ac, the crossing point between rs 
and r-^"^^, is represented by a dashed line. The dotted-dashed line is the outer ergosphere rt • 
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a=0.1M 




a/M (Q/M) 



Stable(L_) Unstable (-L+) 



Q/M (a/M) 



FIG. 20: Black hole case: Last stable circular orbits risco(L-) (thick black curve), risco(—L-) (thick black dashed curve) and risco(—L-\-) 
(thick black dotted curve) as functions of Q/M, for increasing values of the source spin. The angular momentum of an orbiting test particle 
is shown in each region. Black regions are forbidden. In the region r < risco > flsco)^ circular orbits are unstable (stable). The curves 
L = —L— exist in a region located in r < r_. Right: Stability regions at fixed a/M(Q/M) as functions of Q/M(a/M). 
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IV. NAKED SINGULARITIES 
A. Orbits with zero angular momentum 



In this section, we start the study of the motion around a KN naked singularity. We will focus first on the "orbits" 
that are solutions of the following equations 



V\r,Q,a,L) = 0, L = 0. 



(61) 



Note that real solutions can exist only in the case of naked singularities. This particular type of circular motion can 
occur in many axially symmetric spacetimes, and has been investigated in detail for neutral and charged test particles 
in the RN and Kerr spacetimes in [6]-[8]. In the RN spacetime, this particular orbit was interpreted as corresponding 
to a particle located at rest dit r = = Q'^ /M with respect to static observers located at infinity. The radius 
= Q^/M coincides with the value of the classical radius of an electric charge, which is usually obtained by using a 
completely different approach (see also [22 -24 ). In general, the motion defined by the conditions 61, is governed only 



by the gravitational component of the effective potential (and the electromagnetic interaction in the case of charged 
particles), with a zero centrifugal force. Therefore, if a configuration with L = is possible, it should be the result 
of a balance of forces. Thus, this fact can be interpreted in terms of repulsive gravity effects. It is then important to 
note that these limiting conditions can exist only in the case of naked singularities, even for a = or Q = 0. This 
means that the "repulsive component" of the effective mass in the case of a Kerr naked singularity is guaranteed by 



the so-called "spin-charge" of the source. The presence of both charges in the KN spacetime, in Eq. 14, makes this 
effect to become more complex due to the combination of these two components. For more details about the repulsive 
gravitation effects in naked singularities, see also [25 . 

The solutions of Eq. 61 are listed in Table VII and plotted in Fig. 21 , from which it is possible to see that for Q > 




FIG. 21: The radius r/M of circular orbits in Kerr-Newman geometry with zero angular momentum (V' = 0, V = E/ ji and L = 0) 
is plotted as a function of the charge-to-mass ratio Q/M and intrinsic angular momentum a/M oi the naked singularity. The surfaces 
r = r+ and r = r—, delimited by the outer and inner horizon respectively are also plotted. The region r < r+ is shaded. 



there always exists at least one "circular orbit" radius characterized by the conditions given in Eq. 61 for a > a^, and 
in the case of a RN naked singularity; thus, when a = and Q > M, it is r = = Q'^/M. For (5 = 0, "circular 
orbits" defined by Eq. pllexist only in a closed interval a/M G [1, 3a/3/4]. A numerical analysis shows that the orbital 
regions in the plane {a~r) are bounded by the intrinsic spins ai,a2 whose exact form can be found in Appendix [Cj 

This is a feature of the Kerr naked 
le situation is much more complex. 



61 



Moreover, for Q/M < a/5/4 there are several radii satisfying the conditions 
singularity, Q = 0, as well as of the KN singularity. However, in the KN case t 
and the possibility appears to have three radii (ri,r2,^3), but only in the range a < a2. In Figs. 22 we show some 
specific examples of these cases. This analysis is also useful to show that eventually a minimum appears, at which 
these peculiar conditions are satisfied, that corresponds to the interesting case of particle at rest {L = 0) on a stability 
point. 
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Case: Q = (Kerr metric) 


< Q < 0.344263M 


0.344263M < Q < a/5M/4 


a/M Radius 


a/M Radius 


a/M Radius 


M<a<3V3M/4 (f2,f3) 
a = 3V3M/4 f2 


as < a < a2 (ri,f2,f3) 
a = a2 (ri,f2) 
a > a2 ri 


as < a < ai r\ 
a = ai (ri.rs) 
ai < a < a2 (^,^2,^3) 
a > a2 ri 


Case: a/5M/4 <Q <M 


Q> M 




a/M Radius 


a/M Radius 


a > as fi 


a = OVM 
a > fi 



TABLE VII: Description of the circular orbits for a test particle in a Kerr-Newman naked singularity with L = 0. 




B. Circular motion and stability problem 

In this section, we focus on general circular motion with L ^ 0. First, it turns out to be convenient to analyze 
separately the case Q > M, in which the spacetime metric has a singularity not covered by the horizon for any value 
of the spin a/M, and the case Q < M, in which the source spin must comply with the naked singularity constraint 
a > as. In order to simplify the discussion it is useful to analyze first the case Q > >/3M; in fact from Figs. 



23 



and Figs. 32 - which will be discussed below -we expect for Q > ^/3M a forbidden region only for r < r*, and only 
three orbital regions namely (we use a classification similar to that introduced in Sec. |III[ ): region I characterized by 
L = L_, region II with L = {L_, —L-^)^ and region III with L = —L±. A fourth orbital region IV with L = — L_ 
appears for M < Q < a/3M. The characterization of circular motion for these values of charge-to-mass ratios is 
considerably more complex and, therefore, requires a different approach. 

Thus for Q > ^/3M and < a < a^ circular orbits exist within the interval < r < f. The explicit value of 



f depends on the ratios a/M and Q/M (see Fig. 23). It can be shown that in this interval, the orbital angular 
momentum is L = —L±. On the other hand, for r > f there exist circular orbits with both L = L_ and L = — 
For sufficiently large values of Q/M, all these circular orbits are stable. The last stable circular orbit (risco) is, 
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therefore, located at r = r*, where the particle angular momentum and the energy are 



respectively. Note that this was exactly the case for a = 0, in which for Q/M > >/5/2 all the orbits with r > 
were stable, while for lower values of Q/M we had a closed region of instability. The energy is always positive and 
independent of the value of the intrinsic angular momentum of the gravitational source; it approaches the value 



^ fi for large values of Q (see Figs. 23). 



For a > circular orbits exist at r = with angular momentum and energy as given in Eq.(62). For the 
remaining radial distances there are several possibilities. In fact, in the interval r^c < r < f, the orbits have L = —L±, 
in r < r < rs the orbital angular momentum can be either L = — or L = L_, in ra < r < r4 it is L = L_, and, 
finally, for r > there are circular orbits with L = L_ and L = — L+. The situation for Q > V^M is summarized in 
Table |VIII[ In Appendix [Bj the particular case Q = 2M is studied as an application of the results obtained for this 

Case: Q > VSM 



Region 


Angular momentum 


a > aT 


]QVM,r] 


-Li 




{L-, -L+) 


]^3,r4] 


L- 


]r4, oc[ 


{L-, -L+) 


< a < ttT 


]QVM,f] 


-L± 


]f,oo[ 


{L-,-L+) 



TABLE VIII: Description of the circular orbits for a test particle in a KN naked singularity with Q > VSM. 



case. 



The case requires a different analysis because there are some important details in the dynamics around 



naked singularities that must be emphasized. In Figs. 23, (for Q > M) and Figs. 24 (for Q < \/3/2M), we plot the 



orbital radius in terms of the intrinsic angular momentum a/M for different values of Q/M. For a better presentation, 
we have "rotated" the figures, plotting the inverse of the function a = a{r), and using the definition in Eq. (45). We 
want to emphasize that in the plane (r, a) , these curves coincide with f so that the correspondence between the results 
in Table VIII and Figs. [23] does not give rise to confusion. Thus the highlighted regions in Figs. 23 are the orbital 
regions with fixed Q/M for different values of the spin. In all the cases, the region r < is forbidden. For Q > ^/3M 
the situation is rather simple: As one can see in Figs. [23j and it is also shown with details in Table VIII[ there are 
three allowed orbital regions which we have denoted as follows: I characterized by L = L_, II with L = (I/-, — ^+), 
and III with L = —L±. Using this classification, we can interpret the results in Table VIII as follows: The region 
III has a finite, bounded spatial extension that reaches a maximum Tmax at a given value of the spin amax^ and then 
decreases as the spin of the source increases. The values {vmax, ctmax) depend in general on the source charge-to-mass 
ratio Q/M. The size of this region increases with increasing Q/M. The regions I and II, instead, have an infinite 
extension. 

The situation for M < Q < \/3M is much more complex and in addition to the three orbital regions, which 
characterize the motion in Q > V^M, there is now a fourth orbital region IV with L = — L_. The size of the orbital 
region IV increases as Q/M decreases. Instead, the region III decreases with decreasing Q/M until it disappears in 
the limit Q = M (see Fig. [23|). The orbital region IV increases until it reaches a maximum and then decreases. The 



figure highlights the presence of a further type of naked singularity for < (9/8)M^; in fact, for = (9/8)M^, an 
additional banned orbital zone appears for a < and r < r+. This zone increases as Q/M increases, and for fixed 
values of a/M and Q/M it is localized in a finite spatial region. Notice that at r > r+ there are orbits for spin a < ai 



IX 



we include the different cases. 



with (I/-, — ^+), that is, there exists a region of type II. In Table 

Finally for Q < M, and a > ag (condition of existence of a naked singularity with Q < M, Eq. (21), the analysis 
is schematically summarized in Figs. [24] Some examples are discussed in Appendix. [C] We note that the region IV, 
with L = — L_, has a very complex structure and, in contrast to the case. M < Q < V^M, it not only possesses a 
maximum at Vmax^ but also a minimum Vmin that appears as the charge decreases. This situation is illustrated for 
Q = 1/2M in the central plot of Fig. 24 We see that, for a suitable range of values of the source angular momentum. 
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Q>M 



Q=M 




FIG. 23: Left upper: The radius f/M is plotted as a function of the momentum a/M. The dots represent the points {aT,r(aT}- For 
< a < ttT circular orbits exist in the region r* < r < f with L = —L±, for r > r with L = L_ and L = — 1/+ (see text). Center: The 
orbital angular momentum L(r^) / (fiM) and (inset) the energy E(r^)/ ji are plotted as functions of Q/M within the range [1, 10]. Upper 
right and bottom: Plots of the intrinsic KN angular momentum a_/M (gray curve), a+/M (black curve) and aijM (dashed gray curve) 
as functions of r/M, for increasing values of the charge-to-mass ratio Q/M. The angular momentum of the orbiting test particle is shown 
in each region. Black regions are forbidden: circular motion is not possible there. Dotted lines represent the the radii r*, and ry 
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TABLE IX: Description of the circular orbits for a test particle in a Kerr-Newman spacetime. At r 



singularity case for Q < M, it is always a > ag 
and Figs. [32] 



And tti < tts in r G]r^ [5 whereas ai 



as for and r 



Ty it is a+ = d—. In the naked 
M. See also Appendix [a| 



as the radial distance increases there appear two disconnected regions of type IV. This situation becomes extreme 
with decreasing values of Q/M; for instance, for Q = 1/lOM the two regions are completely separated. The outer 
region, farthest from the singularity, has a finite extension and the inner region, closer to the singularity, has an 
infinite extension. 

The radius at which the region II, with starts, increases as the source charge-to-mass ratio Q/M 

decreases. This can be interpreted as a consequence of the fact that the effects due to the spin become increasingly 
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Q-O.IM Q=0.5M Q=(V3"/2)M 




FIG. 24: Plots of the intrinsic KN angular momentum a—/M (gray curve), a+/M (black curve) and ai/M (dashed gray curve) as 
functions of r/M, for increasing values of the charge-to-mass ratio Q/M = 0.1,0.5,0.866025. The region a > as is explored: a naked 
singularity occurs for Q < M and a > as- The angular momentum of the orbiting test particle is shown in each region. Black region 
r < r* is forbidden; no circular motion is possible there. The dotted-dashed line is the outer ergosphere rt ■ 



predominant as the charge decreases, affecting the properties of the counter-rotating orbits with L = — L+. As a 
result, circular motion is no more allowed at very close to the central source where the effects of the gravitational 
stress due to the singularity are expected to be very strong. This effect is corroborated by the fact that a similar 
behavior is found when the charge remains fixed and the spin of the source is increased. 

It is important to notice that once the two parameters (a/M^Q/M) are fixed within in a suitable interval, the 
regions where circular orbits are allowed become disconnected. In particular, for Q < M and orbits sufficiently close 
to the singularity, the regions affected by this phenomenon are of the types IV and I. 

As a result of this analysis, we can identify four types of naked singularities according to the charge-to-mass ratio: 

NS-I : For Q/M > VS with orbital regions of the type I, II and III 

NS-II : For 3/2^2 <Q/M <V3 with orbital regions of the type I, II, III, IV and one forbidden region 
NS-III : For 1< g/M < 3/(2^^2 with orbital regions of the type I, II, III, IV and two forbidden regions 
NS-IV : For Q/M < 1 with orbital regions of the type I, II and IV 

Notice that in naked singularities of the type NS-IV, it is necessary to consider the position of the ergosphere 
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region [r7,r+], because, as shown in Figs. 

The results of studying the problem of stability are represented in Figs. [25] The open curves shown in Figs. [2Q| for 
black holes, find a smooth continuation here in the naked singularity case, as illustrated in Figsj25] As a result, for 
fixed values of Q/M, a finite region of instability appears that splits the plane into two regions of stability. The plots 
show a disconnected structure that resembles the structures found previously in the case of RN [6 and Kerr [7] naked 
singularities. 

Several features appear in the case of naked singularities. Firstly, there is an new region of exclusively counter- 
rotating orbits (type III) that is practically always stable. Also, in general, the radii of the orbits risco intersect now 
at a point, thereby creating a complex structure which depends on the source spin and charge. The region I, which 
is closer to the singularity (say at r < 5M), presents the most complex structure so that by varying a/M, we can 
have regions completely filled with either stable or unstable orbits or also partially filled with stable and unstable 
orbits. In this regard, for small enough spin (say approximately a < lOM), the region I splits into two disconnected 
regions that become larger as a/M decrease, until they become entirely separated by a type IV region {L = as 



illustrated in Figs. 33 This region appears to be always unstable and is bounded by the curve risco{L-). Then, the 
region of instability, which increases with increasing a/M, will eventually cover at a = lOM the single type I region, 
which therefore results in an almost complete island of instability. For sufficiently large values of a/M, the radius 
risco{L-) > r+. 

In conclusion, the existence, extent and stability properties of the region I is mainly determined by the spin of the 



source. We clarify the situation in a schematic picture of the stability region in Figs. 25 (bottom-right). 
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FIG. 25: The charge-to-mass ratio Q^/M^ Q^/M and as functions of r/M^ for increasing values of the source spin. The angular 
momentum of the orbiting test particle is denoted in each region. Black regions are forbidden. The radius of the last stable circular 
motion risco is plotted: risco(L-) (thick black curve), risco{—L-) (thick black dashed curve) and risco{—L^) (thick black dotted curve). 
In the region r < risco (^ > f^lsco)^ circular orbits are unstable (stable). Bottom right: An approximate schematic illustration of the 
stability regions (see text). A naked singularity occurs for a > M and Q > Qs] then, Qs = 0.994987M for a = O.IM, Qs = 0.707107M 
for a = 1/V2M, Qs = 0.6M for a = 0.8M, Qs = 0.43589M for a = 0.9M, Qs = 0.141067M for a = 0, 99M, and = for a = M. The 
dotted-dashed curve is the ergosphere boundary re- 
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V. BLACK HOLES VS NAKED SINGULARITIES 



This section summarizes the main results concerning the circular motion and stability properties obtained for the 



black hole case in Sec. Ill and for the naked singularity case in Sec. IV We show a benchmarking and discuss the 
key differences between these two cases. We start by considering the extension of the stability analysis performed 
in Sec. |TlT] and illustrated in Figs. [19] to include the naked singularity region, as shown in Figs. [27[ The left plot 
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FIG. 26: Left: The radii r^^^^ = risco(—L-\-), ^^^^o ^ ^isco(L-), and = i^lscoi—L-) of the last stable circular orbit in the Kerr 

spacetime as functions of the source angular momentum a The horizon r+ and the radii (dotted line), rj" (dashed line), rs (black 
curve), and r4 (gray curve) are also plotted. The black region (r < r+) is forbidden. The value ac, at the crossing point between rs and 
^Zsco' represented by a dashed line. The gray banned regions correspond to the naked singularity case (NS), the white ones to the 
black hole case (BH). Right: The last stable circular orbit radius risco in the Reissner-Nordstrom spacetime for a black hole Q < M and 
a naked singularity Q > M with different charge-to-mass ratios Q/M. 



shows the case Q = (Kerr spacetime) discussed in Sec. Ill and extensively studied in [T. For all values of a/M, 
the last stable circular orbit risco{—L-\-) exists, that is, for a naked singularity it is the continuous extension of the 
radius for the black hole case. On the other hand, the radius the last stable circular orbit in the black hole case 
with momentum L = L_, does not extend to the values a > M. Then, in accordance with [7 , two additional radii 
determining the stability of circular motion appear: the radius risco 

for 1 < a/M < 3^3/4 and risco{L-) for 
a/M > 3V^/4. The introduction of a source charge changes qualitatively the stability pattern. Firstly, the behavior 
of the two radii with domain only in the naked singularity region drastically depends on the value of Q/M. This can 
be seen by comparing the left and right plots of Figs. [27| A similar behavior characterized the case a = 0. Indee d, as 

for sufficiently small charge source 1 < Q/M < \/9/S, 
< Y^5/2 the curve risco changes direction so that for 



III 



and in Fig. 26 

< Q/M 



shown in the naked singularity case in Sec. 
there appears only one risco ^ while in the range 1 
Q/M > Y^5/2 all allowed orbits are stable. 

The stability analysis of the naked singularity case must be carried out in detail for orbital regions near the 
singularity. We will do this by using the analysis performed in Sec. |IV| for a naked singularity, by extending the 
results of Figs. 25 to the black hole case, and by matching Figs. 20 and Figs. 25 We will also use all the quantities 



introduced in Sec. |IIB| Thus, the analysis will be carried out in parallel in two alternative ways: (1) We reinterpret 
the analysis in Figs. |26|by studying the function r = r(a): by fixing a particular value of the charge-to-mass ratio, we 
vary the source spin in order to cover both the black hole and the naked singularity regions. However, it turns out to 
be convenient to use the inverse function a = a(r) of Figs. [26j The obtained regions are the same as in Figs. 26, but 
with a suitable enlargement of some regions close to the singularity. This is the generalization of the work performed 
in [7 to include a charge-to-mass ratio. This analysis turns out to be very useful to understand the main difference 
between the allowed regions for the circular motion around the two sources. (2) In the second analysis, we fix the 
source spin a/M and vary the charge Q/M to study the orbits r = r{Q). This is a generalization of the case a 7^ of 
the left plot of Fig. 26 for the RN spacetime. This allows us to determine how the rotation of the source affects the 
results of the plot. This is the generalization of the work performed in [6 to take into account the spin-to-mass ratio. 
It turns out also in this case to be convenient to use the inverse function Q = Q{r) in order to explore the details 
of the orbital regions. This comparative analysis is very important to study in detail the stability properties of the 
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FIG. 27: The radii 



o{-L+), r- 



flsco{—L—) of the last stable circular orbit in a KN spacetime in 



terms of the source spin a/M, with different values of charge-to-mass ratio. The horizon r+ and the radii (dotted line), rj" (dashed 
line), rs (black curve), r4 (gray curve) are also plotted. The black region (r < r+) is forbidden. Left: For the orbits with L = — L_ there 
is a stable and an unstable region bounded by the radius tisco{~L—)\ for r < tisco{~L—) the orbits are unstable and for r > tisco{~L—) 
are stable. . 



circular motion around black holes and naked singularities. 



A. Comparative analysis of the circular motion 



We now summarize the results of comparing the properties of circular motion in the entire KN spacetime. It is 
convenient to confront our results with those obtained for the limiting cases a = and Q = 0, which were presented 
in [6-8 , and are here schematically represented in Figs. 26 It was found that black holes and naked singularities are 



characterized by connected and disconnected regions of stability, respectively. The analysis of the generalization that 
includes the spin and charge of the source simultaneously is shown in Figs. 28 -31 for a varying charge-to-mass ratio 



for selected and increasing values of the spin a/M . Figures 32 show the behavior of the circular orbits in terms of the 
spin a/M, for selected and increasing values of the charge Q/M. 

We first consider the analysis of Figs. 28 jSo} The regions of existence of circular orbits with their orbital moments 
are marked in white and gray, while the black regions are those in which the motion is not possible. In Fig. 28, we 
consider the particular case a = O.IM and explore the circular orbit properties for both the black hole and naked 
singularity fields. Notice that from a physical point of view, in the case of black holes, we are interested in the 
motion outside the outer horizon only, r > r+, which in the plots corresponds to the range [r+, +oo[U[0, a^] U [0, Qs]- 
Nevertheless, to understand better the behavior of the relevant radii and for the sake of completeness, we include 
some plots of spatial regions located very closed to the origin. 

In Figs. [29] and [30j we consider the cases a = 0.8M and a = M, respectively. We see that the location of the 
regions where circular motion is allowed depend drastically on the value of the ratios a/M and Q/M. 

We summarize the analysis in Fig. [3l] where we show the behavior of the orbital radius for increasing values of the 
source angular momentum. As pointed out in Sec.[lV| it is possible to distinguish four different regions where circular 
motion is allowed with different orbital angular momenta: The region I with L = L_, II with L G — III 
associated to L = —L± and, finally, the region IV with L = — L_. The region r < (or equivalently Q > Q*) is 
forbidden for both black holes and naked singularities. A second forbidden region appears for configurations with 
a < M and Q < Q~ . This region becomes smaller as a/M approaches the value a = M, where it disappears. We 
note the peculiar value a/M = 1/V2 which determines the interval 1/a/2 < a/M < 1, where the two forbidden 
regions are completely disconnected and separated on the left by a region with L = — L_ and on the right by a region 
with L = L_. In other words, for a < 1.299M there exist two closed regions with L = L_ which are disconnected, 
becoming connected in a particular region for larger values of a/M. The region with L = —L± exist within the interval 
Qc < Q < Q^' The open range with L = — and L = L- begins in the region r > Vx- Furthermore, the region with 
L = — L_ becomes smaller as a/M increases to become practically nil as the area of the region L = —L±. The region 
with orbits L_ becomes increasingly wider as the intrinsic angular momentum and the charge of the source increase. 
We conclude that the structure of the orbital regions around a KN source is very complex in the intervals a < M and 
r <rM2- 
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FIG. 28: The case: a = O.IM. Properties of circular orbits in a KN spacetime, including black holes and naked singularities. The 
charge-to-mass ratio Q^/M, Q^jM and Qc-2 are plotted as functions of r/M. The value of the angular momentum of the orbiting test 
particle is shown in each region. Black regions are forbidden; no circular motion is possible there. Dotted lines represent the radii fo, ''^M2 5 
and r*. The top left plot shows the range r/M G [0, 4]. The remaining plots magnify the intervals where non-trivial behaviors are observed. 
The solid (dashed) white curve is the outer (inner) horizon r+ (r_). We denoted with the lowest value of f(Q = 0). Th horizontal 
thick black line represents Q = Qs- A naked singularity occurs for a > M and Q > Qs, i-e., Qs = 0.994987M in this case. The black hole 
(BH) and naked singularity (NS) regions are denoted explicitly. The dotted-dashed curve represents the ergosphere boundaries r^. 
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FIG. 29: The case: a = 0.8M. Properties of circular orbits in a KN spacetime, including black holes and naked singularities. The 
charge-to-mass ratio Q^/M, Q^/M and Qc2 are plotted as functions of r/M. The angular momentum of the orbiting test particle is 
shown in each region. Black regions are forbidden; no circular motion is possible there. Dotted lines represent the radii fo, rM2, t^ms, t*, 
rx- The top left plot shows the range a/M G [0,5]. The remaining plots magnify the intervals where non-trivial behaviors are observed. 
The solid (dashed) white curve is the outer (inner) horizon r+ (r_). We denoted with the lowest value of r{Q = 0). The horizontal 
thick black line represents Q = Qs- A naked singularity occurs for a > M and Q > Qs, i.e., Qs = 0.6M in this case. The black hole (BH) 
and naked singularity (NS) regions are denoted explicitly. The dotted-dashed curve represents the ergosphere boundaries r^. 
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FIG. 30: The case: a = M. Properties of circular orbits in a KN spacetime, including black holes and naked singularities. The charge- 
to-mass ratio Q^/M, Q^/M and Qc2 are plotted as functions of r/M. The dotted-dashed curve represents the ergosphere boundaries . 
The angular momentum of the orbiting test particle is shown in each region. Black regions are forbidden; no circular motion is possible 
there. Dotted lines represent the radii fo, r*, r^. The left top plot shows the range a/M G [0, 6]. The remaining plots magnify the intervals 
where non-trivial behaviors are observed. The solid (dashed) white curve is the outer (inner) horizon r+ (t— ). We denoted with the 
lowest value of f(Q = 0). A naked singularity occurs for a > M and Q > Qs, i-e., Qs = M in this case. 
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FIG. 31: The charge-to-mass ratio as a function of the radial coordinate for increasing values of the source angular momentum. The 
ratios /M, Q^jM and Qc-i are also plotted. The angular momentum of the orbiting test particle is shown in each region. Black regions 
are forbidden; no circular motion is possible there. Dotted lines represent the radii fo, r*, and r^- We denoted with the lowest value of 
f(Q = 0). The horizontal thick black line represents Q = Qs- A naked singularity occurs for a > M and Q > Qs, i-e., Qs = 0.994987M for 
a = O.IM, Qs = 0.707107M for a = 1/V2M, Qs = 0.6M for a = 0.8M, Qs = 0.43589M for a = 0.9M, Qs = 0.141067M for a = 0, 99M, 
and Qs = for a = M. The black hole (BH) and naked singularity (NS) regions are denoted explicitly. The dotted-dashed curve 
represents the ergosphere boundaries . 
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As seen in Sec. |III[ an alternative way to study the neutral test particle circular motion in a KN spacetime is 
to explore the orbital regions, for different values of the ratio a/M, as functions of the radial coordinate. This is 
particularly suitable in order to compare the results with those of the Kerr spacetime depicted in the left plot of 



Fig. 26 In Fig. [32j we show the behavior of the circular orbits for different values of the charge-to-mass ratio of the 
source. In this case, fixing the source charge and varying its intrinsic angular momentum, we plot the regions where 
circular motion occurs and point out the value of the orbital angular momentum associated with each region. The 
boundaries of the regions are defined by the intrinsic spin a± and a^. For Q < M there are essentially three types of 
orbital regions characterized respectively by L = — L_, L = L_, and L = (i^-, — ^+). We can see that there are two 
forbidden regions, the first determined by the conditions r < and a < for Q^/M^ < 9/8, and the second one by 
r < for Q^/M^ > 9/8. The region within the interval L+) increases as Q/M increases, while the regions 

with L = — L_ and L = L_ decrease as Q/M increases and, finally, they coincide at r = r*. Interestingly, the region 
with L = — L_ disappears for Q'^/M'^ > 3. 

One can see that the orbital region (L = L_) diminishes as the parameters Q/M and a/M decrease, and that 
it extends from a rather small distance from the central source (light gray region). For Q/M small enough, say 
Q <C M/2, this region disappears and, as the charge increases, it moves towards areas with higher spin and smaller 
radii. For Q > M a new orbital region appears, characterized by two different types of counter-rotating orbits with 
L = —L±] this region becomes increasingly larger as the ratio Q > M increases. Finally, we note that the regions 
with L = —L± and L = —L- are bounded and extend up to a maximum radius and a maximum angular momentum 
of the source. 

The curve a = ai divides the plot into two regions, the left black one, in which no circular motion is allowed, and 
the right one, in which orbital motion occurs for a > with L = L_ and for a < ai with L = (i^-, — ^+). As the 
source spin increases, the left forbidden region first grows to a maximum value of the orbital radius and then begins 
to decrease. On the other hand, the region with L = (i^-, — ^+) (right gray region) covers larger orbital radii as the 
spin increases, and moves towards the central source as the charge-to- mass ratio decreases. 
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FIG. 32: The intrinsic KN angular momentum a-/M (gray curve), a+/M (black curve) and ai/M (dashed gray curve) as 
functions of r/M, for increasing values of the charge-to-mass ratio Q/M . The angular momentum of the orbiting test particle 
is shown in each region. Black regions are forbidden; no circular motion is possible there. Dotted lines represent the radii r*, 
and r^. The dotted-dashed white curve is the outer horizon r+. The dotted-dashed black curve represents the ergosphere 
boundaries r^. The thick black horizontal line is a = a^; a naked singularity occurs for Q < M and a > a^. The black 
hole (BH) and the naked singularity (NS) regions are explicitly denoted. In this case, as = 0.994987M for Q = O.IM, 
as = 0.866025M for Q = 0.5M, and a, = 0.5 for Q = V3/2M. 
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B. Discussion of the stability problem 



In this last section, we will study the influence of the intrinsic angular momentum and charge of the source on 
the stability properties of test particles, moving along circular orbits around the central source. In Sec. |IIIB| and 
Sec. |IVBl we studied the case of a black hole and a naked singularity, respectively. In Fig. 33, we show a review of the 
results when both sources are analyzed simultaneously. First, we note the presence of last stable circular orbit radii 
risco for all possible orbits with angular momentum — L+, — L_ and L_. The location of risco depends on the value 
of the intrinsic angular moment. The larger is the rotation, the further from the origin is risco- This is physically 
reasonable since the rotation is known to generate its own gravitational field in general relativity. On the other hand, 
for a fixed value of a/M, as the charge-to-mass ratio increases, the value of risco first increases until it reaches a 
maximum value and then it decreases. We have, in general, that risco{L-) > risco{—L-\-) > risco{—L-)^ but this 
situation becomes significantly more complex with increasing spin and decreasing source charge. 

To formulate the main results in a plausible manner, let us imagine an idealized accretion disk made of test particles 
which are moving along circular orbits on the equatorial plane of a KN spacetime. It then follows that the structure 
of the accretion disk will be determined by the locations of the radii of the last stable circular orbits. Consider, 



for instance, in Fig 33 the case a = 0.99M. The radius risco{L_) determines the minimum radius of the idealized 
accretion disk which in the case of a black hole could be extended, in principle, to infinity, because the stability zone 
has no upper bound. The situation is different in the case of a naked singularity. The disconnected structure of the 
stability regions is especially clear for a > M. However, for a = 0.99M, for example, the radius r^sco(— ^+) clearly 
goes into the NS zone creating a stable (left) orbital region with L = — L+, an unstable region limited by risco{—L-^), 
and a stable region (right) for r > risco{—L-\-). This situation does not occur in the black hole case. Notice a second 
unstable region of type L_ separated by a stable zone (right) from a —L_ region. 
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FIG. 33: The charge-to-mass ratios Q^/M, Q^jM and cis functions of r/M^ for increasing values of the source spin. The angular 
momentum of the orbiting test particles is shown in each region. Black regions are forbidden. The radius of the last stable circular 
motion risco is plotted: risco(L-) (thick black curve), risco{—L-) (thick black dashed curve) and risco{—L+) (thick black dotted curve). 
In the region r < risco (^ > t^Isco)-/ circular orbits are unstable (stable). Bottom right: An approximate schematic illustration of the 
stability regions (see text). The horizontal thick black line represents Q = Qs- A naked singularity occurs for a > M and Q > Qs, 
i.e., Qs = 0.994987M for a = O.IM, Qs = 0.707107M for a = 1/v^M, Qs = 0.6M for a = 0.8M, Qs = 0.43589M for a = 0.9M, 
Qs = 0.141067M for a = 0, 99M, and Qs = for a = M. The black hole (BH) and naked singularity (NS) regions are denoted explicitly. 
The dotted-dashed curve represents the ergosphere boundaries r^. 
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VI. CONCLUSIONS 



In this work, we performed a detailed analysis of the behavior of test particles moving along circular orbits on the 
equatorial plane of the KN spacetime. We use the method based on the investigation of the effective potential from 
which one can derive the location and the stability properties of the circular orbits. We also analyzed the energy and 
angular momentum of all possible orbits, and used this information to identify different orbital regions which allow 
us to find out the differences between black holes and naked singularities. 

We have identified four orbital regions, namely: region I characterized by L = L_, region II with L = {L_^ —L-^), 
region III with L = — and a fourth orbital region IV with L = — L_. The region I contains only co-rotating 
particles whereas regions III and IV are for counter-rotating particles. It is interesting that region II contains both 
co-rotating and counter-rotating particles with different orbital angular momentum. 

We found that naked singularities can be divided into four classes that depend on the value of the charge-to-mass 
ratio: NS-I for {Q/M > V3) with orbital regions of the type I-II-III, NS-II for (3/(2^2) < Q/M < V3) with 
orbital regions of the type I-II-III- IV and one forbidden region, NS-III for (1 < Q/M < 3/ (2^/2)) with orbital 
regions of the type I-II-III- IV and two forbidden regions and NS-IV for {Q/M < 1) with orbital regions of the type 
I-II-IV. Figure 34 illustrates this classification also for the case of black holes. 
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FIG. 34: Classification of KN compact object according to the charge-to-mass ratio Q/M and the spin a/M of the source. Black holes 
can have regions of type I and II and only naked singularities can accommodate regions of all four types. The type NS-III (with *) has 
two forbidden orbital regions while NS-II has only one forbidden orbital region. 



We found that from the point of view of the circular motion, there exists always a forbidden area (r < r*) which is 
completely inaccessible for test particles. In the case of black holes, the radius is situated inside the outer horizon. 
In the case of naked singularities, it can be considered as a surface that "covers" the singularity itself and prevents it 
from being studied with orbiting test particles. The value of the radius = Q^/M does not depend on the source 
angular momentum so that it is a property of the gravitational field generated only by the electric charge only . It is 
interesting that the expression for coincides with the classical radius of an electric charge Q with mass M which is 
obtained by applying a completely different approach. In the limiting case of a RN naked singularity, this radius has 
a very particular significance because it corresponds to particles with vanishing angular momentum (L = 0) [6 . 

Another interesting result is that only naked singularities allow orbiting test particles of the type III and IV, i.e., 
counter-rotating particles with L = —L±. Since black holes can have counter-rotating particles only with L = — L+, 
it follows that the presence of counter-rotating particles with L = — L_ would necessarily imply the existence of KN 
naked singularities. 

In the case of black holes, we have shown that there exists only one stability region for circular orbits whereas 
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naked singularities present two stability regions separated by one zone of instability. This peculiar structure suggests 
the possibility of distinguishing black holes from naked singularities by investigating the structure of accretion disks 
around the central source. Indeed, if we imagine an idealized accretion disk made of test particles which are moving 
along circular orbits on the equatorial plane of a KN spacetime, then in the case of a black hole the accretion disk is 
continuous whereas in the case of a naked singularity it is discontinuous. It would be interesting to find out if this 
difference will remain in the case of more realistic accretion disks j26ti32] . and if it could lead to observable effects. 
This question was also recently analyzed in [33 . 

During the past few years, much effort has been made to understand the question about the formation of naked 
singularities [34ll36] , and about the properties of solutions that describe the gravitational field of such hypothetical 
compact objects [37^44] . The present work contributes to this discussion in the sense that it shows the possibility of 
identifying KN naked singularities by investigating the properties of their accretion disks, the KN spacetime. 

However, there are indications that the rigid structure of the KN solution is not adequate to correctly describe 
naked singularities. In fact, several studies show that KN singularities are highly unstable configurations [45, |46] 
that, therefore, cannot exist in Nature. It follows that it is necessary to consider more general solutions that could 
describe the gravitational field of compact objects. The first possibility is to consider solutions with quadrupole 
and higher moments [47 . Since already in the case of a static field with only quadrupole moment, there are many 
possible solutions [48H5^ and most of them are quite difficult to handle, we propose to use the recently proposed 
reinterpretation of the Zipoy-Voorhees [55^-^ solution which from the mathematical point of view is the simplest 
generalization of the Schwarzschild black hole solution with a quadrupole moment. We expect to investigate this 
problem in the near future. 
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Appendix A: Limiting cases 



To study circular motion in the gravitational field described by the KN spacetime, in Sec. II B we introduced the 
spin parameters a^, and a± which determine the boundaries of the regions where circular motion is allowed for 
both black holes and naked singularities. When these parameters coincide, the extended regions allow only a limited 
number of orbits which depend on the charge-to-mass ratio of the central source. In this Appendix, we present the 
characteristics of these limiting orbits. 

First, orbits with a+ = a_ are characterized by the radii 
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1 

3M 
1 

3M 



+ -2^ sin 



2ysm 




arccos 



y 



M V 



3M 
1 

6M 



y 



1 /3 

2<3^ + 3\/3 (8 - 13Q2 + i6Q4)^ 
+ 3v^VQ6 (8 - 13g2 + IQQ 



1/3 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 

(A6) 



with y = ^/M^ - +Q^, and fj = 8M^ - UW^Q^ + ISM^Q^ + gQ^. The behavior of these radii is illustrated 

in Fig. [35 
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Moreover, for the limiting case = we obtain the following radii 
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Finally, the condition = a_ leads to the orbits characterized by the radii 
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In Fig. p5| we show the behavior of the orbits in terms of the ratio Q/M . In Fig. 35 left, we consider the first set 
of orbits {r7/2, r*, Ttt, f^i, f2y2, T'yg}. These radii exist for all values of Q/M] they decrease when Q/M decreases and 
diverge for as Q/M increases. This means that the ranges where these radii have boundaries become larger as Q/M 



increases. The situation is completely different for the radii of the second set, {vp 



plotted in Fig. 
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right. These orbits determine a closed region of the space r x Q. The radii r^, where = ax, and where 
ai = a_, are defined for a closed range of values of Q/M and are bounded. The distance between these two radii can 
become null and has a maximum as a function of the charge. 




FIG. 35: Left: Plot of r*/2, r* 



and ry^ as functions of Q/M. Right: Plot of , , and Vp as functions of Q/M. 



Appendix B: Naked Singularities with Q > ^/SM 

Particular cases for the naked singularity with Q > ^/SM NS-IV are sketched and discussed in Figs. [36l[37l[38l[39l 
[40|[4lJ|42j which show the effective potential and the orbital angular momentum for different values of the spin a/M. 
We explore circular motion of neutral particles in a KN naked singularity with Q = 2M, where arj^/M = 46/3V^, 
analyzing the orbits in the regions < a/M < 46/3^3, a/M > 46/3 V^, and a/M = 46/3 V^. 
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a/M a/M 

FIG. 36: Stability regions for a KN naked singularity with charge Q = 2M (left) and Q = AM (right). The radius f (thin black line), rs 
(thick gray line) and r4 (thick black line) are plotted as functions of the naked singularity spin a/M. The gray region r < r^^ = Q'^ / M 
is forbidden. The spin ax is denoted by a thick gray dashed line. The inset plot is a zoom. White regions are stability zones. In the 
light-gray regions the orbits with angular momentum — 1/+ are unstable. 
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FIG. 37: Naked singularity with charge Q = 2M and intrinsic angular momentum a = AM . Left: The energy V{—L—) of (stable) circular 
orbits with angular momentum L = — L_ is plotted in the region r* < r < r. The radius r* (dashed line) and f (dotted line) are explicitly 
represented. Center: The energy V{L—) of (stable) circular orbits with angular momentum L_ is plotted in the region r > r. Right: The 
energy V{—L-^) of (stable) circular orbits with angular momentum L = — 1/+ is plotted in the region r > r*. The radii r* (dashed line) 
and r (dotted line) are shown explicitly. The inset plots are a zoom in the region r* < r < r. 
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FIG. 38: Naked singularity with charge Q = 2M and intrinsic angular momentum a = 12M. Left: The energy V{—L-) of (stable) 
circular orbits with angular momentum L = — L_ is plotted in the region r* < r < f. The radii r* (dashed line) and r (dotted line) are 
shown explicitly. Center: The energy V(L—) of (stable) circular orbits with angular momentum L_ is plotted in the region r > r. The 
radiii r4(dotted-dashed line) and rs (gray line) are plotted. The inset plot is a zoom in the region r < r < rs- Right: The energy y(— L+) 
of circular orbits with angular momentum L = —Lj^ is plotted in the region r > r*. In the black region ra < r < r4 no orbits with angular 
momentum L = — 1/+ exist. The inset plot is a zoom in the region r* < r < rs. 
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FIG. 39: Circular orbits in a KN naked singularity with Q = 2M and different values of a/M. The picture plots the locus of critical 
points of the effective potential V/ ji with different values of the momentum parameter a for a < Q (left plot) and a > Q (center and right 
plot). The (orbital) angular momentum, L/(Mfi), of circular orbits is plotted as a function of the orbital radius r/M. In the central plot 
the regions r/M G [4,10] and L/^Mji) G [—70,0] are explored. The right plot focuses the details on the zoomed region r/M G [4,4.30] 
and L/{Mii) G [-5,1]. 
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FIG. 40: Circular orbits in a KN naked singularity {Q = 2M) are explored. Left: The radius f/M is plotted as a function of the 
momentum a/M . The effective potential, the orbital angular momentum and the energies are plotted as functions of r/M. 
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FIG. 41: Circular orbits in a KN singularity with Q = 2M and a = lOM are explored. The energy £'(zbL-t), the angular momentum 
±L±, and the effective potential V/ ji of particles in circular orbits are plotted as functions of the orbital radius r/M. 
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FIG. 42: Circular orbits in a KN naked sing ularity with Q = 2M and a = 46/(3v^)M are explored. The energy E{^L±) the angular 
momentum zbL-t and the effective potential V/^ of particles in circular orbits are plotted as functions of the orbital radius r/M. 
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Appendix C: Naked singularities with Q < M 

In Sec. IV A, we use the expressions for the spin ai = a_ and a2 = a+ with 
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in the study of solutions of the equations of motion with L = 0, where 

^5 = 3 (27 + 128Q^ - 48Q^f - 192Q2 (144 - 81Q^ + 92Q^ + 32Q^) - 

2592 21/3Q4 (9216 - 10368Q2 + 15529Q'^ + 2496Q^ + 5232Q^ + 2560Q^^) 

^3 

3 (19683 - 93312Q2 _ 337392Q4 + 2522240Q^ - 3866496Q^ + 780288Q^o + 4096Qi 

§22/3^3 + 1024^4 



22/3 



^3 



(CI) 



(C2) 



L4 



1024^3 



1728 2^/^Q'^ (9216 - 10368Q^ + 15529Q'^ + 2496Q^ + 5232Q^ + 2560Q 



(729 + 32Q2 (-72 + 593Q^ - 568Q'^ + 8Q^)) ^3] ; 



/3 



(C3) 
(C4) 



i2 = 93312Q^ (884736 - 1492992Q^ + 2656080Q'^ - 819693Q^ + 1183756Q^ + 426790Q^°- 
+ 924864Q^^ + 33280Q^'^ + 65536Q^^) ; 



(C5) 



Li = 34828517376Q^° (lO - 37Q^ + 16Q'^)^ (8 - 13Q^ + 16Q'^) (783 + 31ig^ + 128Q'^) 



(C6) 



To facilitate the presentation we use in the above formulas the normalized quantities Q ^ Q/M and r r/M. The 



behavior of the spin parameters is depicted in Fig. 43 




FIG. 43: Circular orbits in Kerr-Newman naked singularity for Q < M are explored. We plot the parameters ai, 0,2 and ao as functions 
of the charge Q/M. 
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T hese p arameters were used in Table VII to determine the existence regions for solutions with L = 0. As discussed in 



it is important to consider in particular the ranges Q/M g]0, 0.344263] with (ai, a^), Q/M g]0. 344263, >/5/4[ 



with (ai,a2,as), and Q/M G [a/5/4, 1] with only. 
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To study explici t exam ples, we considered the case Q = 0.8M in Figs. |44|45[ Q = 0.5M in Figs. [46| |47|48[ [49) and 
Q = 0.2M in Figs. 5Q|51 We use throughout the expressions ri = r~_^ r2 



+ '• 



The plots show the behavior of the effective potential for negative and positive values of the orbital angular 
momentum. These plots confirm the presence of a minimum for the effective potential also with L = 0; thus, there 
are critical points of the effective potential associated with stable configurations for the test particle orbits. We notice 
also the presence of minima with negative energy values. 



Q^O.SM 



Q=Q.8M 




FIG. 44: The KN naked singularity with Q = 0.8M and a > \/3/2M ~ 6M. Left: The radii r, ri, r2, r4 as functions of the intrinsic 
angular momentum a/M. Left-inset: Plot of the (orbital) angular momentum, L/(Mfj,) as a function of the orbital radius r/M, for 
different values of a. Right: The radii of last stable circular orbit risco(L-)^ (thick black curves) and risco(—L-\-) (thick black curve) are 
plotted as functions of a/M. The radii f, ri, r4 are also plotted. The inset plot focuses on the region r/M G (0.6, 2) and a/M G (2, 11). 
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FIG. 45: Circular orbits in the KN geometry with a = 1.5M and Q = 0.8M are explored. The effective potential is plotted for different 
values of the angular momentum L. The radii r, ri, and r4 are plotted (dashed lines). The critical points of the effective potential are 
represented by points. 
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Q=0.5M Q=0.5M Q=0.5M 




r/M r/M a/M 

FIG. 46: Circular orbits in a KN naked singularity with Q = 0.5M and a > a/3/2M ^ 0.866M are explored. Left and central plots: 
Plots of the locus of critical points of the effective potential V/fi with different values of the momentum parameter a for a > 0.866M. The 
(orbital) angular momentum, L/{Mii)^ of circular orbits is plotted as function of the orbital radius r/M. Right: The radii r4^/M^ r^/M^ 
r/M^ fi/M^ r2/M f^/M are plotted as functions of a/M. 
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FIG. 47: Circular orbits in a KN naked singularity with Q = 0.5M and a = I.IM are explored. Left: The effective potential for 
L = -A.ZMjji, L = -3M/i, L = -M/i, L = -Mji, L = Mfi. For L = -4.3M/i there is a minimum at r = 0.261M with energy 
E = -2.227/^, a second minimum at r = 11.314M with E = 0.965/i and a maximum at r = 7.356M with E = 0.966/i. For L = -3M/i 
there is a minimum at r = 0.264 with E = — 1.519/i. For L = —Mji there is a minimum located at r = 0.287M with E = — 0.283/i. For 
L = 0, the minimum is at r = 1.013M with E = 0.326/i, for L = M/i there is a minimum at r = 1.808M with E = 0.676/i. Notice the 
presence of the minima with negative values of the potential. Center and Right: The energies E(^L±) and angular momentum of circular 
orbits as functions of r/M. 
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FIG. 48: Left: Circular orbits in a KN naked singularity with Q = 0.5M and a = 2M are explored. The effective potential for L = —5Mfi, 
L = —2Mii and L = Mfi is plotted. For L = —bM/j. there are two minima, one minimum is at r = 0.254M with energy E = — 1.341/i 
and the second one at r = 17.01M with E = 0.975/i, there is a maximum at r = 8.353M where E = 0.99/i. For L = — 2M/i there is a 
minimum at r = 0.257 where E = — 0.32/i. For L = M/j, a minimum is located at r = 0.304M where E = 0.843/i, a maximum at r = 0.568 
with E = 0.848/i. There is an additional minimum at r = 2.446M where E = 0.781/i Notice the presence of the minima with negative 
values of the potential. Center and Right: The energies E(^L±) and angular momentum of circular orbits as functions of r/M. 
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FIG. 49: Circular orbits in a KN with Q = 0.5M and a = 1.3M are explored. The effective potential is plotted for L = —AAMfi^ 
L = — 2M/i, L = Mfi and L = 0. Notice the presence of the minima with negative values of the potential. Center and Right: The energies 
E{^L±) and angular momentum of circular orbits as functions of r/M. 




r/M 

FIG. 50: Left: Circular orbits in a Kerr-Newman naked singularity with Q = 0.2M, and a >^ 0.979796M) are explored. Left: Plots 
the (orbital) angular momentum, L/{Mii) as function of the orbital radius r/M, for different values of a. Center: A Kerr-Newman naked 
singularity with Q = 0.2M, and a > ag. The radius fi, f2, rs, ri, and r4 are plotted as functions of the intrinsic angular momentum a/M. 
Right: A Kerr-Newman naked singularity with Q = 0.2M, and a > 0.6M. Right: The radii of the last stable circular orbits risco{L-)^ ( 
black curves) and nsco(~-^+) (thick gray curve) risco(—L-) (gray curve), are plotted as functions of a/M. The raadius r2,^3, ^4 are also 
plotted. The inset plot focuses on the region r/M E (0, 1) and a/M E (1, 1.5). 




FIG. 51: Circular orbits in the Kerr-Newman geometry with a = M and Q = 0.2M are explored. The effective potential is plotted for 
different values of the angular momentum L. The radii fi, f2, fs, ri, and r4 are plotted (dashed lines). The critical points of the effective 
potential are represented by points. 
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